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CONTINUUM  THEORIES  FOR  CYLINDRICAL 


WAVE  PROPAGATION  IN  A LAMINATED  MEDIUM 
by 

Yakc  Benveniste 

Department  of  Solid  Mechanics  Materials  & Structures 
School  of  Engineering 
Tel -Aviv  University 
Tel-Aviv  Israel. 

Abstract 

In  this  report j interacting  continuum  theories  are  developed  for 
cylindrical  wave  propagation  in  a bi-iemlr.ated,-,  linearly  elastic,  composite 
medium  We  consider  a cylindrical  circular  infinite  cavity  perpendicular 
to  the  layering  of  the  ccmpcsite . The  laminated  medium  is  subjected  at 
the  cavity  wall  to  a non  axisymmecric  shear  and  normal  loading  uniform  in 
the  direction  of  the  cavity  axis,.  Two  microstructure  theories  are 
developed:  the  first  simulates  the  effectively  two-dimensional  plane 

motion  cf  the  medium  which  is  due  to  a normal  ana  shear  loading  in  the 

circumferential  direction  the  second  models  the  effectively  anti-plane 

*'  * '-T 

motion  due  to  a shear  leading  in  the  direction  of  the  cavity’s  axis,  \ 
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1.  Introduction 

The  investigation  of  wave  propagation  in  a composite  medium  is  usually  carried 
out  by  constructing  a model  of  the  composite  which  describes  in  some  sense  its  gross 
mechanical  response.  Such  an  approach  is  unavoidable  since  an  "exact"  description  of 
the  mechanical  response  oi  the  composite  is  almost  impossiDle  due  to  the  complex 
geometry  and  the  multiple  reflections  that  waves  undergo  at  tne  interfaces  of  its 
constituents-  The  simplest  model  or  a composite  medium  is  the  sc-calxed  "etfectlve 
modulus  theory"  which  replaces  the  composite  dv  an  equivalent  homogeneous  anisotropic 
material  with  equivalent  elastic  constants  and  descnoes  the  gross  mechanical 
behavior  in  terms  of  the  averages  of  displacements,  stresses  and  strains  over  the 
representative  elements.  Ihe  effective  modulus  theory,  however,  is  not  able  to 
account  for  a very  important  phenomenon  in  the  propagation  of  waves  in  composite** 
that  is,  dispersion  and  attenuation. 


C 


A more  realistic  description  of  the  mechanical  benavior  of  composite  media  can 
be  attained  by  taking  into  consideration  in  some  sense  its  microstructure.  Several 
such  models  exist  now  in  the  iiteracure  and  recent  review  papers  by  Achenbach  [1] 
and  Moon  [2]  cover  the  majority  of  the  research  works  carried  out  in  this  area. 

One  approach  to  model  a composite  by  taxing  into  account  the  microstructure  i*  to 
use  mixture  theories  as  models  of  composites  as  suggested  oy  Lempriere  [3].  These 
are  interacting  continuum  theories  vhere  the  constituents  are  superimposed  in  space 
bat  allowed  tc  undergo  individual  deformations-  Tne  basic  difficulty  in  these  theories 
is  the  analytical  specification  of  the  interaction  terms  which  arise  m the  formulation 
For  the  specific  case  of  longitudinal  wave-guide  type  propagation  in  laminated 
composites,  Hegemier  et.  al.  [4]  use  an  asymptocic  method  whose  first  order  term  can 


Numbers  in  the  brackets  designate  references  at  the  end  of  tne  paper. 


be  cast  in  :ne  form  oi  a binary  mixture  theory,  and  give  a rational  construction  of 
both  mixture  interaction  and  constitutive  relations.  By  comparison  with  exact  and 
approximate  phase  velocity  data,  they  show  that  the  developed  theory  gives  good  first 
mode  agreement  for  wavelengths  greater  than  the  typical  composite  micro-dimension, 
furthermore  good  correlation  with  experimental  data  in  transient  wave  propagation  is 
reported  in  £4]  „ The  validity  of  this  theory  m transient  wave  propagation  was  also 
checked  by  Aboudi  by  comparing  results  with  those  obtained  by  solving  che  complete 
dynamic  equations  of  elasticity  ra . In  Nayfeh  and  Gurtman  extended  che  theory 
to  study  shear  waves  in  laminated  wave  guides.  Again,  very  satisfactory  1 -suits  ate 
obtained  by  comparing  the  dispersion  relations  with  exact  phase  velocity  spectra  for 
the  ’..cwest  mode  of  propagation.  Ihe  theories  developed  in  (V)  and  are  effectively 

one-dimen sionai  in  rhe  sense  that  the  mixture  variables  are  dependent  on  one  space 
variable  only. 

In  this  paper  continuum  microstructure  theories  of  the  type  given  in  [4]  and  [7] 
are  developed  for  che  case  of  two-dimensional  wave  propagation  from  a cylindrical  cavity 
in  a bi-iaminaced  medium.  We  consider  a cylindrical  circular  infinite  cavity  perpen- 
dicular to  the  layering  of  the  composite  and  define  a cylindrical  coordinate  system 
(r.  0,  z)  with  z being  along  the  axis  of  the  cavity  (See  Fig.  1).  Continuum  mixture 
theories  will  be  developed  for  two  leading  situations  at  the  cavity  wall-  a)  Situations 
in  which  the  loading  is  tangential  in  the  ©-direction,  and/cc  normal  to  the  cavity^ 
with  dependence  on  0 , but  uniform  in  the  z direction,  and  b)  Situations  in  which 
the  loading  is  tangential  In  the  z-direction  with  dependence  on  0 but  again  uniform 
in  the  z-direction.  Clearly,  a specific  example  of  a ©-dependent  and  z- independent 
loading  is  an  infinite  line  load  in  che  z-direction.  Symmetry  considerations  shew  that 

case  (a)  yields  symmetric  u ^ and  u ^ and  anti- symmetric  u ^displacements  with 

t y z 

(si)  (W) 

respect  to  the  midplanes  of  each  layer,  whereas  (b)  yields  anti- symmetric  1^.  & and 

1 I 

ft  i 
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(a)  * 

symmetric  uz  displacements  again  with  respect  to  the  midplanes.  It  is 

(ot) 

clear  that  when  the  displacements  are  averaged  over  a layer  no  net  u displace- 

z 

ments  are  obtained  in  case  (a)  and  no  net  u^  and  displacements  in  case 

r o 

(b) . Therefore  we  will  characterize  the  motion  in  (a)  by  the  name  quasi-plane 

motion  and  that  in  (b)  by  the  name  quasi-anti-plane  motion.  Note  that  both 

motions  are  two-dimensional.  For  the  specific  situations  in  which  the  loading  at 

the  cavity  wall  is  axi-symmetric , three  cases  are  distinguished:  (1)  Axi-symmetric 

(a) 

normal  loading  at  the  cavity  wall;  in  this  case  there  are  no  u^  displacements 
(a)  (a) 

and  the  u^  and  u^  displacements  are  respectively  symmetric  and  anti- 


symmetric with  respect  to  the  midplanes.  Since,  on  the  average  the  u 


(a) 


displacements  vanish,  this  motion  can  be  characterized  as  quasi-radial  motion. 

(2)  Axi-symmetric  tangential  loading  in  the  6-direction  at  the  cavity  wall; 

in  this  case  there  are  no  u*"0^  displacements  and  the  ujj0^  and  u^  are 

r 6 z 

respectively  symmetric  and  anti-symmetric  with  respect  to  the  midplanes.  The 
motion  is  characterized  as  quasi-rotatory  motion.  Note  that  cases  (1)  and  (2) 
are  special  cases  of  the  quasi-plane  motion.  (3)  Axi-symmetric  tangential  loading 
in  the  z-direction;  the  motion  caused  by  this  kind  of  loading  is  a special  case 
of  the  quasi-anti-plane  motion  which  does  net  have  any  6-dependence. 

In  the  second  section  of  this  paper,  the  equations  of  motion  and  constitutive 
equations  for  each  constituent  will  be  given.  In  the  third  and  fourth  sections, 
continuum  microstructure  theories  will  be  developed  for  the  quasi-plane  and 
quasi-anti-plane  motion  respectively. 


Here  and  in  the  sequel  the  superscript  or  the  subscript  will  take  the  values 
1 or  2 and  will  indicate  that  the  quantities  belong  to  either  one  of  the 
constituents. 
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2 • Basic  Equations 

Consider  a periodic  array  of  two  alternating  isotropic  linearly  elastic  layers 
of  widths  2h^  and  21^  respectively.  Let  us  also  consider  an  infinite  circular 
cylindrical  cavity  perpendicular  to  the  layering.  Define  a cylindrical  coordinate 
system  (r,  9,  z)  such  that  the  z coordinate  is  along  the  axis  of  the  cylindrical 
cavity.  Furthermore,  let  zq  be  a local  coordinate  measured  from  the  midplane  of 
each  layer  (See  Fig,  1). 


The  equations  of  motion  and  the  constitutive  equations  in  cylindrical  coordinates 
for  every  individual  layer,  are  given  by  : 


+ h “if + h -if  * i “{?>*  j - p.  fo  ■ia) 
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(2.1) 
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where 


.(a) 


3_  (a) 
ar  ur 


ul  + [ ( Tq  u 


<L  <«) 

ae  ue 


uia))/r  ] 


i_u(a) 

dza  2 


(2.10) 


and  and  u^^  with  i and  j standing  (or  r,  6,  z are  the  stresses  and 

displacements,  p is  the  density  of  each  constituent,  X and  u are  the  Lame 

oc  a a 

parameters  and  t denotes  the  time. 

3.  Quasi-plane  Motion 

In  this  section,  we  will  consider  situations  lr.  which  the  uva,>  displacements 

z 

are  anti-symmetrlc  and  the  u^  and  u^"1 1 displacements  symmetric  with  respect  to 
the  midplanes  of  the  layers.  The  general  case  of  a c-dependent  ncrmax; and  shear 
loading  in  the  0-direction  at  the  cavity  wall  wixl  be  treated  and  a continuum 
mixture  theory  will  be  developed  which  wili  modex  the  composite  in  this  two- 
dimensional  quasi-plane  motion  Ac  the  end  of  the  section,  the  special  case  of 
axi-symmetric  wave  propagation  will  be  briefly  discussed. 

To  develop  the  micrcstructure  model  the  equations  of  motion  (2.1)  and  (2,2) 
will  be  averaged  over  the  z^  coordinates  and  alter  defining  "partial  stresses" 
and  "partial  densities"  they  will  be  cast  in  a standard  binary  mixture  form.  A similar 
averaging  process  will  be  carried  out  for  the  constitutive  equations  (2,4),  (2,5)  and 
(2.7)  and  approximate  constitutive  relations  connecting  the  partial  stresses  to  the 
average  constituent  displacements  wixl  be  obtained 


We  define  the  average  quantity 

h 


ijJ<a)  (r  ,0,t)  - (1/h  ) 
a 


a 


^aj(r,0,z3,t)dza 


(3,1) 


A 


and  the  partial  stress  and  densities  as: 


.o» 

a r9 


h'2 


n s(a) 
a 90 


na  " ha^h  » with  h " h1  + h2>  (3.2 

Note  that  the  symmetiy  of  the  and  Ug3^  and  the  anti-symaetry  of  u^°^ 

displacements  imply  symmetry  with  respect  to  the  midplanes  of  the  a^a\  0^a\ 

rr  99  zz 

stresses  and  anti-symmetry  of  and  . We  now  take  the  average  of 

equations  (2.1)  and  (2.2), and  employing  the  stated  anti-symmetry  and  continuity 
conditions  of  the  0^  and  0^  stresses  at  the  interfaces  we  obtain  the  following 
mixture  equations  of  motion: 

h «/*>  Iff  ♦ K»SP>-  '$>>M  - .wi  ;<“) . 


qaQ(r,0,t) 


(3.4) 


"fr  °r0P^  + (1/r)  haQ0P)+  (2ar0P)/r)  ' P uQ(a)=  qaP(rt0,t)  (3.5) 

® t* 


where 


k 
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[Q(r  ,0,c)  ] (h)  = O^(r,0,h9,t)  = c/^(r,0,-h  ,t)=  -o^(r,0,h  ,t) 

is*  s.  is*  x rz  x 


(3.6) 


[P(r,6,t)](h)  = O^(r,0,h2)t)  = (^(r.O.-h^t)  = -a^J  (r.O.h^t)  (3.7) 


and 


qi  = 1 


q2  = -1 


We  will  call  the  partial  stresses 

rr 


P , 0 P and  the  average  displacements 
tv 

-(a)  -(a)  0y 

, Uq  by  the  term  mixture  variables.  The  next  step  is  to  develop  approximate 


constitutive  equations  relating  the  partial  stresses  to  the  average  constituent  dis- 
placements. Furthermore,  expressions  relating  P and  Q to  the  mixture  variables 
will  be  obtained.  These  expressions  will  be  based  on  the  assumption  that  for  a 


characteristic  wave  length  L , and  composite  microdimension  h = h^  + h^  , the  ratio 


£ = h/L  is  a small  parameter  and  terms  of  order  e can  be  neglected  (see  [4]). 


Averaging  equations  (3.4  - 3.7)  according  to  (3.1)  and  using  the  anti-symmetry 


(a) 


and  continuity  properties  of  u we  obtain: 


[(0(a)/A  ) - (E  /X  ) u(a) 
rr  a a a 9r  r 


(1/r) ( §q-  u£a)  - u£a))]na  = -q  S(r,0,t)  (3.8) 


;(a) 


,9  - (a)  -(a),  3 -(a). 


Iloee  'V  - 


(3.9) 


- 'W  fe“t“)'(VV(1/rK  le  "ea>  " ■ -VV.e.t) 


'’r“P)  ■Vc.<(1/rt(l?"r“)-:;9a)>+l7;;e“>' 


(3.10) 

(3.11) 


In  order  to  abbreviate  the  equations,  this  notation  will  be  used  in  the  sequel  too. 
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[Q(r,e,c)](h)  = a^)(r,0,h2,t)  = ^(r.e.-^.t)-  -o^)(r,0,hJ  ,t) 


(3.6) 


[P(r,0,t)](h)  = o^)(r,0,h2,O  = (r.e.-hj^.t)  = -a^J  (r.S.h^t)  (3.7) 


ql  = 1 


q2  1 ' 


We  will  call  the  partial  stresses  o 


(ap)  (ap)  (ap) 
rr  ’ Ooc  ’ °r0 


and  the  average  displacements 


-(a)  -(a) 

u^_  , Uq  by  the  term  mixture  variables.  The  next  step  is  to  develop  approximate 

constitutive  equations  relating  the  partial  stresses  to  the  average  constituent  dis- 
placements. Furthermore,  expressions  relating  P and  Q to  the  mixture  variables 
will  be  obtained.  These  expressions  will  be  based  on  the  assumption  that  for  a 

characteristic  wave  length  L , and  composite  microdimension  h = h^  + h2  , the  ratio 

2 

e = h/L  is  a small  parameter  and  terms  of  order  e can  be  neglected  (see  [4]). 

Averaging  equations  (3.4  - 3.7)  according  to  (3.1)  and  using  the  anti-symmetry 

(cx) 

and  continuity  properties  of  u^  we  obtain: 

l(5<“>rt0>  - (E^/y  i-C™  - (l/r)  ( ;‘a>  - Gf,)Jna  - -V(r,8,t>  (3.8) 


[(59?'V  - <VVa/r)(!r  "e<0l)+"'O)> 


(3.9) 


‘<5«'V  - <W  h ^“>-<VE«)a/r)(  h “»a)  + ^a,)K  - 


(ap)  r/1/  . , 3 -(a)  -(a).  . 3 -(a). 

°r0  “ Va  [(1/r)(  36  ur  " u0  > + u0  ] 


(3.10) 

(3.11) 


In  order  to  abbreviate  the  equations,  this  notation  will  be  used  in  the  sequel  too. 
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where 


E 

a 


A 

a 


+ 2u 


(3.12) 


and 


(h) [S(r ,0,t)] 


(r,0,h2,t) 


(r,0,-h1,t) 


u(1> 


(r.U.h^.t) . 


13.13) 


We  will  now  obtain  an  approximate  expression  for  S , which,  when  substituted 

in  (3.8)  and  (3.9)  will  give  together  with  (3.11)  the  desired  mixture  constitutive 

(cx)  (cx)  (rx) 

equations.  For  this  purpose  we  expand  u , ' and  u in  turns  of 

utilizing  their  symmetry  properties.  Using  then  these  expressions  in  (2,6)  an  expansion 
( cO 

for  a is  obtained  (See  [4]).  This  is  of  the  form: 
zz 

a(c°  = 5(a) + 0(e2)]  (3.14) 

zz  zz 

Using  the  continuity  of  across  the  interfaces,  up  to  the  order  of  approximation 

we  are  concerned  with  in  this  paper,  we  can  write: 


(3-15) 


Equation  (3.15)  together  with  (3,10)  furnishes  an  expression  for  S 


S(r,9,t)  = (A1/E)[  |^-  u^1}  + (l/r)(u<i}  + |g-  u^)]  - 


- (V»i|r“r2>  T a/r)(G‘2)+|j;<2>,  , 


(3.16) 


where 


E • (E2/n2)  + (E.^/^) 


Substituting  (3^16)  in  (3^8)  and  ('3,9)  we  obtain 
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(3.19) 


a00P^  = c22^1/r;(  90  u0  ^ T u J2)>1  + 1*  l s7  + (i/i)(  fgSj0*  i?>i  1 


b -(2) 

f d,  -r-  u 
2 9t  t 


»J»> 


where 


c m [n  E - (a2  /E)  ] 
aa  a a a 


Ca$  » XaA„/E 

da  * <Vc  - '»> 


a -2D 


with  0,8  • 1,  2 and  a 4 8 


■ ■ . : 


- 1C  - 


Equations  (3  io)  and  (3  17-3  20,  axe  the  ar-air^a  mixture  ..nditux.ve  squat  or. 

It  remains  now  to  aerive  expressions  relating  ^ c*na  r to  one  mixture  vaxj.aDx=s 
Lee  us  first  obtain  tne  expression  ror  Q Aga^n  cr,c  pro.^dure  icuowsa  in  Laj  *Ui 
Dfc  used  here  and  it  wij.1  be  given  oriei^  A^th^ut  abi^w  Mui.i..r ^ equation  ,2  b, 
Dy  , expanding  j.n  powers  of  ana  »ntegr=.tAng  Dy  parte  we  .ou^n: 


ur  (r,o,h  ,t)  - u t ,ti  / 3y  l u '.r,'3,n,L;  - >.1/^  't-j'-sn  ,^/j  — 

1 >*  I X O l Z Ot  C f.  L 2*  \ 


a . u; 


UeJr.g  one  definitions  (3  6)  and  (a  ij;  g*v£s. 


a^(r,b,h ,,t)  - + vn.h,3>i-  — or  Qj-0 

i l r x dr  * 

u^2)(r,9,n2,tj  - u^2>  r (r..n,  3;  L — o - U/n2>  Qj -C 


Subtracting  (3  24)  from  (3  23) , employing  the  .cntmu_ty  condition  „f  a^a;  act  . . 
the  interfaces j ana  using  equation  (3  *6),  a£  get: 


Q(r,0,t)  = (K/h2KuU) 


)Cu  ) T Jf  iAt(X  / E)  uvi^-  r 

x r . x ( *-  l 

▼ .£;  C.U'  - vv./E)  u^2>]) 

s j.  0 x 3 


whe  r e 


K = [(nj/3yt^)  + (n2/3p2)i 


Vo  1.C) 


M = K/3 


and  and  ^ are  differential  operator  r ,=u~h  that 

ii<*>  * lA* 

Zcw  - 57  * * (l/r)» 

(1/r)  fee  * “ (i"2' 


! 


The  expression  for  P can  be  obtained  by  applying  the  same  steps,  this  time  to  equation 
(2.9).  The  result  is: 


P(r,0,t)  = (K/h2)^1*-  G<2))  + [Jf2  + (2/rZ)  {Ml  (A^u^-CX^E)^  ] } 


y»x-(D. 


+ (1/r)  ~ (M[(A1/E)Gj1)-  (a2/E)u^2)]>  (3.30) 

90 

The  mixture  equations  of  motion  (3. 4, 3. 5)  together  with  the  interaction  terms  Q and 

P given  by  (3.25)  and  (3-30),  and  the  mixture  constitutive  equations  (3.11,  3.17-3.20) 

define  completely  the  quasi-plane  motion  of  the  laminated  medium  as  modeled  by  the  present 

mi crostructure  theory.  Note  that  the  equations  are  two-dimensional  in  r and  0 and  the 

tangential  displacement  and  the  radial  displacements  u£°^  are  coupled  to  each  other. 

Substitution  of  (3.11,3-17-3.20,  3.25  and  3.30)  the  equations  cf  motion  (3.4,  3*5) 

provides  four  coupled  partial  differential  equations  for  u^  , , u^  and  u£2^  ; when 

r r * o u 

these  are  prescribed  at  the  cavity  wall,  together  with  initial  conditions,  they  can  be 

determined  throughout  the.  laminated  medium. 

For  the  special  case  of  axi-symmetric  loading  at  tne  cavity  wall,  there  is  no  depend- 

— (ot)  "“(q.) 

ancc  on  0 and  inspection  of  the  field  equations  show  chat  the  u^_  and  u^  displacements 
decouple  each  other  In  cne  case  there  is  a net  radial  motion  ot  the  medium  and  m the 
other,  a net  rotatory  motion;  these  are  called  respectively  the  quasi-radial  and  quasi- 
rotatory motion  of  the  laminated  composite - 

4.  Quasi-anti-plane  motion 

In  the  present  section  we  will  consider  the  motion  oi  the  composite  which  is  caused 

by  a tangential  loading  in  the  z-direction  at  the  cavity  wall,  this  input  being  uniform 

(ot) 

in  the  z-direction  but  dependent  in  0.  The  motion  will  De  such  that  the  ug  displace- 
merits  are  symmetric  with  respect  to  the  midplanes  whereas  the  u^  and  the  displace- 

ments are  anti-symmetric.  A net  motion  in  the  z-direction  will  be  present  which  will  be 
two-dimensional  in  r and  0.  As  in  the  previous  section  we  will  develop  here  a continuum 
microstructure  theory  modeling  the  laminated  composite  in  tnis  two-dimensional  quasi- 
anti-plane  motion. 


mmam 


- I?  - 


ritsc,  it  is  noted  that  symmetric  u^°^  and  anti-symmetric  u^°^  and  give 

rise  to  svmmemc  a,,a^  and  and  anti-symmetric  C^\  c!^  , a(,^/  , with  respect 

9z  rz  rr  60  zz  r9 

tc  tne  mldpianes  Averaging  equation  (2.3)  according  to(3.1)  and  defining  again  the 

partial  stresses 


C(Q?j  . n 5(a) 
rz  a rz 


gi  ^es 


where 


[ lr  + a/r)jc^f;)+  (i/r)  1 0^p)-  °iP)7T^c,)  = ‘»aR(r»0»t) 

dt 


or 


(4.1) 


(4.2) 


[R(r,0,t)]h  - C1 2) (r ,9,h.,t)  = a(l) (r ,0 ,-h. ,t)  = -a(1)(r,0,  h.  ,t)  (4.3) 

z Z Z ZZ  1 ZZ  1 

Averaging  equations  1 2 . 8 ) ana  (2  9)  we  obtain: 


a('ap')  - y n 4—  u(a)  = q y U(r,6,t) 

rz  a a dr  z a 


a6zP;  " Uana(1/r)  fe  ^z^  = Va  T(r*0»t) 


(4.4) 

(4.5) 


where. 


and 


o<«  • n 5?> 
9z  a 9z 


[U(r,e,t)]h  = ura)(r,9,h1,t)  = i/2)  (r  ,0  ,-h2 , t)  = -u<2)  (r,0,h2>t) 
[T(r,0,t)]h  - u^Cr.e.h-.t)  = u(2)  (r,0,-ho,t)  = -ut2)  (r  ,9,h.  ,t) 

“ JL  / l 


(4.6) 

(4.7) 

(4.8) 


Equations  (4,2>  are  the  mixture  equations  of  motion  and  equations  (4-4)  and  (4.5) 
the  mixture  constitutive  equations.  They  contain  the  functions  R(r,0,t),  U(r,9,t)  and 
T(r,0,t)  which  will  have  to  be  determined  in  terms  of  the  mixture  variables.  The 
determination  ct  these  functions  will  be  carried  out  by  following  similar  procedures 
to  those  given  in  [7],  In  the  present  case,  three  coupled  differential  equations  will 
be  obtained  for  R,  U and  I- 
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Multiplying  equations  (2.6)  by  zfll  expanding  u^°^  and  in  powers 

of  z^  and  then  integrating  by  parts  give  (See  [7]): 

yu<°°  (r,0,ha,t)  - G<“>]  = (ha/3){a^)(r,e,hcx,t)  - 

' V h ura)(r*0«ha*t)  + a/r)u^a)(r,e,ha,t)+(l/r)|QU^a)(r,e,ha,t)]} 

(4.9) 

Using  the  definitions  (4.3),  (4.7)  and  (4.8)  results  in: 

u^1)(r,6,h1,t)  - u^1)  = (h1h/3E1){-R-A1^Pu  - ^(1/r)  T>  (4.10) 


u[2)(r,0,h,,t)  - u(2)  = (h0h/3E0){R  + \2$q  II  + X0 (1/r)  T>  (4. 


90 


11) 


f Ql\ 

Subtracting  (4.11)  from  (4.10),  employing  the  continuity  condition  of  u^  'across 
the  interfaces,  we  get: 


ajtf0  U(r,0,t)  + ai(l/r)  — T(r,0,t)  + a2R(r,0,t)  + [ (u^2) -u^ ) /h2]=0 


where 


(4.12) 


al  = [(X1n1/E1)  + (X2n2/E2)]/3 


(4.13) 


&2  ~ Kn^/Ej^)  + (n2/E2)]/3 


(4.14) 


Equation  (4.12)  is  the  first  differential  equation  relating  U,  R and  T to  u ’ 


Substitution  of  (2.4),  (2.5)  and  (2.7)  in  (2.1)  gives: 

a 


£<“)u(a)  + X™  u<a>  +T2-  a(“>  +X 
^3  r 4 0 9z  rz  a 

a 


9r9z  z 

a 


(a)  n 

u =0 


(4.15) 


where  £<a>  end  }£<“> 


are  differential  operators  such  that, 


(t.M> 


■'A'’1*  Vr  ’ t?  * - fr  * 

0T7  dt 


Jtf’*  * ».  * "X<*>  ‘ «“a/t2)  |f  * 


2,  a 


(ot)  (a) 

Multiplying:  equation  (4.15)  by  x , expanding  u and  u~  in  nowers  of 

Ot  r H 

integrating  hv  narts,  we  obtain: 


(h 


t/^)[4ft)^)(r,e,ha,t)  +^0t)  uJa)Cr,0,hu,t)]  + 


3 


(a) , 


+ x«l7 1 ■ “z"'1  + ^,©,h„>t»  -s "']  - 0 


r(o), 

rz 


(a>! 

rz 


Substitution  of  (4.9)  in  (4.18)  furnishes: 


((no/3)  8t^a)fu^)(r,0>ha,t)/h]}  + {(na/3)^^a;(uj“;(r,0,ha,t)/hl} 


'.»(ot)r  (a) 


where 


+ { (X  n /3K  ) -r—  [o^a\r,0,h  ,t)/h]]  + { [a  ^ (r  ,0  ,h  ,t)  - 

a a a 9r  zz  a rz  a rz 


lla)* 


a?* 


2,  32  S2  , 

— ,I,J  -2-  ft 

30  3t 


- »V'2> 

* ‘(Xa  +V  - (Vr2>  fe* 


Introduction  of  the  definitions  (4.3),  (4.7)  and  (4.8)  in  (4.)0)  provides: 


[ (n^/3)  ^5  ' U]  + [(.n1/3)^i)  T]  - I (A^/31^ ) P]  + 
+ { r°ri><r»®»h1.t)  - o^l/h2}  = O 
-I(n2/3)4tj2)0]  - [(n2/3)^2)Tj  + [{k^/3^}  R]  + 


+ {[o^)(r,e.h,,t>  - o^i/h2}  = o 


(*• IS) 

]/h2)=l 

(4.19) 

(4.20) 

(4.22) 


(4.23) 
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which,  when  subtracted  from  each  other  finally  yields: 


J?7  U(r,e,t)  + #8T(r,e,t)  - ax  ~ R(r,9,t)  + t (o<2p)  /n2) -(o^p)  /n^  ] /h2=0 

(4.24) 


where 


-[<pn1/3)^1)i|»  + (n2/3)oJ^2)^] 

° [(n:/3)  + (n2/3)^2)^] 


(4.25) 


(4.26) 


Equation  (4.24)  is  the  second  differential  equation  relating  U,  T and  R. 


It  remains  now  to  derive  the  third  and  last  differential  equation  for  these 
functions.  We  start  first  by  substituting  (2.7)  and  (2.5)  in  (2.2),  which 
provides : 


+ ♦ tu/r»„  u<“),  + ^ - 0 (4.27) 


•10  r 


a 8z  80  z 
a 


9z  6z 
a 


where 


' 2V1/r2)  !e  l|'  + (Xa  + >V<^2  ^ 


(4.28) 


Applying  to  equation  (4.27)  the  same  procedures  which  were  applied  to  (4.15), 
we  obtain  after  considerable  manipulation: 


J?13T(r,6,t)  +^14U(r,e,t)  - (1/r)  |g-  R + [ (a^2p) /n2) -(a^P^)  ] /h2=0 

(4.29) 


with  and  being  differential  operators  such  that 

Jfl3*  - [(t^/3)^^  +(n2/3)^2)^  ] (4.30) 

- [(t^/3)^*  + (n2/3)^2)^  ] 


(4.31) 


- It  - 


z 2 

£u*  • IEa  (4'32> 

dt 


£u*  ■ «\.+  >c?  - <XX)1(1/t>lri9*  + 


4 i«tt  + 3“0 > - ‘xX)1(1/r2)  h * <4-33> 

Equations  (4.2),  (4.4),  (4.5)  together  with  (4.12),  (4.24)  and  (4  29) 

define  the  quasi-anti-plane  motion  of  the  laminated  medium  as  modeled  by  the 

present  continuum  mixture  theory.  Ihey  constitute  nine  coupled  linear  partial 

differential  equations  for  the  nine  functions  u^  , o^01^  , c^0^,  U,  T and  R, 

z oz  rz 

Since  they  are  a system  of  linear  equacions,  they  can  be  manipulated  so  that  they 
provide  two  coupled  partial  differential  equations  for  and  u ^ » when 

these  displacements  are  prescribed  at  the  cavity  wall  together  with  initial  con- 
ditions, they  can  be  determined  throughout  the  laminated  medium. 


For  the  general  case  of  a 9-dependent  loading  at  the  cavity  wail,  the 

mixture  variables  and  are  functions  of  the  space  variables 

z rz  oz 

r and  6,  as  well  as  time.  For  the  special  case  cf  axisymmetric  loading,  we 

. .1 i _ 1 r.(^p)  ^ ' A . . V „ V.  T _ C *-  U j _ , 1 .. 


expect  that  the  partial  stresses 


should  vanish  In  fact,  this  can  easily 


be  seen  first  by  observing  that  with  ~ = 0 we  have: 

do 

•f2  = = ^14  = 0 


(4.34) 


and  the  terms  containing  T in  (4.12)  and  (4.24)  vanish.  Again  with  = 0, 
equation  (4.5)  gives 


°9z  = Vet  T(r»e*C)’ 


(4.35) 


which,  when  substituted  in  equation  (4.29)  yields  a differential  equation  for  T 
which  is  decoupled  from  all  other  variables.  Thus,  with  an  axisymmetric  tangen- 
tial loading  in  the  z-direction,  clearly  we  have  T=0,  which  through  (4.36)  yields 


17  - 


0,(,ap)=  0.  Equation  (4.2),  (4.4),  (4.12)  and  (4.24)  with  ~ = 0 define  now  the 

UZ  du 

axi-svmmetric  quasl-anti-plane  motion  of  the  composite. 

5.  Conclusion 

Wo  have  developed  an  interacting  continuum  theorv  for  cylindrical  wave  propagation 
■> " a laminated  composite.  The  case  of  a O-dependent  loading  at  the  cavity  wall  is 
treated  and  the  theory  is  two-dimensional  in  the  sense  that  the  averaged  motion  of 
each  constituent  is  dependent  on  the  space  coordinates  r and  0^  as  well  as  .time. 

Note  that  the  construction  of  the  theory  is  strongly  dependent  on  the  specific 
symmetry  conditions  which  arise  due  to  a uniform  loading  in  the  direction  of  the  axis 
of  the  cavity.  Tt  can  be  seen  that  similar  svmmetry  conditions  will  be  present  when 
a laminated  half-space,  with  the  surface  being  perpendicular  to  the  layering,  is 
impacted  hv  a tangential  or  normal  infinite  line  load  perpendicular  to  the  interfaces 
of  the  layers.  Thus,  it  is  expected  that  continuum  mixture  theories  of  the  same  type 
can  again  he  developed  for  these  problems,  these  in  fact  being  the  two-dimensional 
generalizations  of  the  theories  developed  in  [4]  and  [7]. 


L 
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1*  ^ circular  cylindrical  iniirtii*  •■•■iiy  iu  a M-laminated  composite  mediuai. 


